Mixed symmetry localized modes and breathers in binary mixtures of Bose-Einstein 

condensates in optical lattices 



H. A. Cruz^'^, V. A. Brazhnyi\ V. V. Konotop^'^'^, G. L. Alfimov'^, and M. Salerno^ 
^Centro de Fisica Teorica e Computacional, Universidade de Lisboa, 
Complexo Interdisciplinar, Avenida Professor Gama Pinto 2, Lisboa 1649-003, Portugal 
^ Departamento de Fisica, Universidade de Lisboa, 
Campo Grande, Ed. G8, Piso 6, Lisboa 1749-016, Portugal 
^ Departamento de Matemdticas, E. T. S. Ingenieros Industrials, 
Universidad de Gastilla-La Mancha, 13071 Giudad Real, Spain 
Moscow Institute of Electronic Engineering, Zelenograd, Moscow, 124498, Russia 
^ Dipartimento di Fisica "E. R. Gaianiello" , Gonsorzio Nazionale 
Interuniversitario per le Scienze Fisiche della Materia (GNISM), 
Universitd di Salerno, 1-84081, Baronissi (SA), Italy 

We study localized modes in binary mixtures of Bose-Einstein condensates embedded in one- 
dimensional optical lattices. We report a diversity of asymmetric modes and investigate their dy- 
namics. We concentrate on the cases where one of the components is dominant, i.e. has much larger 
number of atoms than the other one, and where both components have the numbers of atoms of the 
same order but different symmetries. In the first case we propose a method of systematic obtaining 
the modes, considering the "small" component as bifurcating from the continuum spectrum. A gen- 
eralization of this approach combined with the use of the symmetry of the coupled Gross-Pitaevskii 
equations allows obtaining breather modes, which are also presented. 

PACS numbers: ... 



I. INTRODUCTION 

Optical lattices are known to be an effective tool for 
manipulation of matter waves. The induced spatial peri- 
odicity allows, in particular, for the existence of various 
localized structures, not available in homogeneous sys- 
tems. Such excitations, which are called localized modes 
(also gap solitons depending on the limit they are con- 
sidered in) have been intensively studied during the last 
few years (see e.g. recent reviews [1, 2] and references 
therein) and observed experimentally in a Bose-Einstein 
condensate (BEG) of ^^Rb atoms [3]. While main atten- 
tion was devoted to single-component BEGs, existence 
of localized modes was found to be a generic property 
of multicomponent systems, as well. More specifically, 
emergence of coupled bright solitons, from the modula- 
tional instability of binary mixtures of BEGs in optical 
lattices (OLs), was observed numerically in [4]. In [5], it 
was shown that OLs can be used to modify the interac- 
tions both within and between BEG components, leading 
to the creation of two-component gap solitons. Existence 
and stability of 2D and ID coupled gap solitons in a bi- 
nary BEG with repulsive interactions was analyzed by 
means of the variational approximation and numerically 
in [6]. It was also found that localized modes can exist 
in atomic-molecular condensates |7|, in Bose- Fermi mix- 
tures in OL d, 0] and in three-component spinor conden- 
sates P!q| . 

All the results on binary mixtures of BEGs reported so 
far, concern excitations of the two components character- 
ized by similar scales, i.e. by spatial localization regions, 
and numbers of particles of the different species having 
the same order. Also earlier studies were restricted to 



modes having density distributions in both components 
of the same symmetry, whenever the components belong 
to the same gap (gap solitons belonging to different gaps 
were addressed in Q). 

In the present paper we report a diversity of novel 
states observable in two-component BEGs embedded in 
OLs with different symmetry properties of the density 
distributions with respect to the periodic potential and 
therefore termed below as mixed symmetry modes (or 
also gap solitons of mixed symmetry). The interaction 
between the two components of the condensate is crucial 
for the existence of these modes since they are formed 
from single component modes which have different sta- 
bility properties (one stable and the other unstable) in 
absence of interaction. The interspecies interaction al- 
lows for stabilization of the unstable mode in the pair 
and subsequently for formation of a stable bound state. 
Mixed symmetry modes exist both for attractive and re- 
pulsive inter- and intra-species interactions and repre- 
sent very stable new excitations of two component BEGs 
in OLs. Stabilization phenomena induced by the inter- 
species interaction are also possible for equal symme- 
try states (both symmetric or antisymmetric with re- 
spect to a minimum or a maximum of the potential). 
In this case we find that modes which are both unstable 
in absence of interaction can form stable equal symmetry 
bound states in presence of interspecies interaction. Be- 
sides stationary states of mixed symmetry we show that 
non stationary states of mixed symmetry which give rise 
to breather-like oscillations are also possible. When the 
inter- and the intra-species interactions are all equal, the 
two-component breather modes can be constructed from 
a linear superposition of two mixed symmetry modes 
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with an unbalanced number of atoms in the two compo- 
nents. On the other hand, we show that breather modes 
with a balanced numbers of atoms and with a very reg- 
ular dynamics which persist on a long time scale with- 
out any apparent emission of matter, are also possible. 
Mixed symmetry modes, both of stationary type and of 
breather type, considerably enlarge the number of stable 
localized excitations which one can find in BEG mixtures, 
this showing the richness of these systems with respect 
to single component BECs. 

The paper is organized as follows. In Sec. [Ill we re- 
call the system of coupled Gross-Pitaevskii (GP) equa- 
tions describing the mixture, specify the parameters to 
be used for particular numerical simulations, and give a 
simple physical picture suggesting existence of asymmet- 
ric modes. In Sec. IIIII we discuss on particular exam- 
ples the localized modes of mixed symmetry in the limit 
{Ni ^ N2) and study stability of the modes by direct 
numerical simulations. In Sec. II VI we consider large am- 
plitude localized modes both of equal symmetry and of 
mixed symmetry type. Sec. |Vl is devoted to breathing 
modes, both with balanced and unbalanced number of 
atoms, which are obtained by combining the developed 
approach with the symmetry of the coupled GP equa- 
tions. In the Gonclusion the main outcomes are summa- 
rized. 



II. STATEMENT OF THE PROBLEM 

A. The model 

At the zero temperature and a large number of atoms, 
when quantum fiuctuations can be neglected, a diluted 
binary mixtures of BEGs is well described by the coupled 
GP equations (j = 1,2) [nj 
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■A*,+T/,(r)*, 



M \7^"^-^l*il' + "i.3-il*3-il'j*i- (1) 

Here rrij is the atomic mass of the j^^ component, M = 
mim2/{mi +777.2) is the reduced mass, and 

(^zj — ^ j i are 

the s-wave scattering lengths of the binary interactions. 
The macroscopic wave functions, ^j, are normalized to 
numbers of atoms Nj-. j\^j\^dY = Nj, which do not 
depend on time. The trap potential for the component 
j, lj(r), consists of a superposition of a magnetic trap 
and an OL along the x-direction: 



the respective recoil energy Erj = / {2mj), and U{-) 
is a TT-periodic dimensionless function (the lattice profiles 
are considered to be the same for both components). 

We consider cigar-shaped condensates, i.e. the limit 
A <C 1, at low densities. To specify the scaling we intro- 
duce a small parameter of the problem e = na, where 
a is an average transverse oscillator length defined as 
a = Y^a/(2MaJ) , where uo = y/cJiCc^ • Then the applicabil- 
ity of the theory is determined as smallness of a number 
of particles on the scale of the scattering length, i.e. as 
Y^n|aii| ~ A>:a = e <C 1, where n = N\f\/a is the linear 
atomic density and N = Ni -\- N2 is the total number of 
atoms. This is the case when the multiple-scale analysis 
can be applied to reduce ([T]) to a system of coupled ID 
GP equations (see e.g. [T,^^)! By introducing slow dimen- 
sionless independent variables X = ex /a and T = e^uot 
and dimensionless order parameter <I>j = C^(l), through 
the expansion (j = 1,2) 



\/2|aii| 

with the normalized transverse eigenfunctions 
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we obtain a system of coupled ID GP equations 



.d<!>j 



nij dX"^ 



+ (Xil*,f +Xl*3-,f)*i. (4) 

Here Voj = V^M/rrij and the nonlinear coefficients have 
the following form 



and X 



Looking for stationary solutions ^je "^^^ , we 



(mi +m2)VcJicj2 ai2 
(micji + m2Cc;2) |aii| ' 

-ifijT 



obtain the system 



rrij dX'^ 



(Xil$il' + xl*3-,f)$i (5) 



with /ii and /i2 being the chemical potentials of the first 
and second components, respectively. 



B. The main goal and the terminology 



rrijujj 



{ri + yx^) + VoERjU{Kx). 



(2) 



Here ooj is the transverse linear oscillator frequency of 
the component j, A is the aspect ratio (it is assumed to 
be equal for both components), d is the lattice constant, 
hz = n/d, Vq is the amplitude of OL measured in terms of 



The goal of this paper is to present some families of lo- 
calized modes and breathers of Eqs. ([5]) which are typical 
for multicomponent BECs and which posses interesting 
symmetry properties. These solutions have been found 
numerically by means of different methods. The typical 
difficulty which arise in this problem is to get a priori 
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idea on possible shapes of localized modes which can be 
used for numerical iterative algorithm. To overcome this 
difficulty we have used two alternative approaches. The 
first involves a search of bifurcations for specific solutions 
of ([5]) using shooting method (see [12]) as the main tool. 
The second is the self-consistent approach which has been 
used before for applications to single-component BECs 
(see [13j). The self consistent method, indeed, is very 
convenient for exploring the band structure in presence 
of nonlinearity and to address to specific modes for which 
good initial guesses are available. Basically, the shooting 
method allows one to make a systematic exploration of 
the whole families of solutions in terms of a single (for sin- 
gle component BECs) parameter (shooting parameter). 
Similar exploration of the case of coupled GP equations 
is quite complicated, since in this case one has to scan a 
semi-infinite domain of two shooting parameters. In the 
present paper we use the shooting method to investigate 
localized modes of the mixed symmetry type in the limit 
A^i N2 (see Sec. Hn]) and resort to self consistent ap- 
proach to investigate BEC mixtures with balanced num- 
ber of atoms (see Sec. lIVp . When one of the components 
is dominant, i.e. it has much larger number of atoms than 
the second component, the last one can be considered as 
bifurcating from the linear spectrum. Such modes will 
be referred to as unbalanced or strongly asymmetric in 
number of atoms. In this way after constructing the lo- 
calized solution for the first component in the absence of 
the second one, the mode can be constructed by means 
of some iterative procedure. Combining this approach 
with the symmetry of the system we will also be able to 
construct (numerically) exact breathing localized modes. 

Another regime we will deal with corresponds to the 
situation when the both components are large amplitude 
solutions having numbers of particles of the same order, 
Ni ~ N2. In this case we will show existence of not 
only simplest fundamental modes with equal symmetry 
but also modes and breathers with mixed symmetry. We 
remark that mixed symmetry modes were also found in 
atomic- molecular condensates In this case, however, 
the mean-field equations had an intrinsic asymmetry due 
to the different properties of atoms and molecules. 

First of all it is desirable to establish a convenient ter- 
minology for the identification of modes. To this end we 
observe that multiplying ([5]) by d^j / dX and integrating 
with respect to X one arrives at the necessary condition 
for existence of a stationary localized mode (<I>i, ^2)'- 



i*ir + vb2i$2i 



dU{X) 
dX 



dX = 0. 



(6) 



For the physical interpretation of this condition, we no- 
tice that a solution of ^ minimizes the Hamiltonian 
H = Hi+H2 + Hi2+ H,{0), where 




dX 




(7) 



H12 = xj |^iPl^2p<^^, and we have defined 

H,{C) = J {Voil^if ^Vo2\^2f )U{X - C)dX. (8) 

Considering (^i(X - C), ^2(X - C)) with C < 1, which 
for a stable solution must lead to increase of the energy, 
we immediately obtain 



dC 



and 



d^H.iC) 



dC 



>0. 



(9) 



The first of the obtained constrains is exactly Eq. ([6]). 

Let us now require the potential to have a definite 
symmetry, i.e. to be either even or odd, and concen- 
trate on the solutions possessing a specific symmetry, i.e. 
^j{-X) = ±^j{X). If the potential is odd, U{-X) = 
— /7(X), such solutions strongly localized about X = 
cannot satisfy (|6|). This allows us to conjecture nonexis- 
tence of stable solitons localized about the point X = 
for odd potentials. In particular, in the case of cos-like 
potential 



U{X) = 2cos(2X), 



(10) 



explored below in the present paper: localized modes 
with a defined symmetry can be excited only in the vicin- 
ity of one of the points X^ = Trn where n is an integer. 

Considering solutions localized about X = for an 
even potential U{X) = U{—X) one can identify four 
types of the modes. They can be classified in terms of the 
symmetry properties acquired with respect to the peri- 
odic potential, considered earlier for a atomic-molecular 
BEC and for a single-component BEC [l^. 

The diversity of modes makes it is desirable to estab- 
lish a convenient terminology for their identification. To 
this end we classify states according to their symmetry 
properties with respect to the minima of the periodic po- 
tential around which are localized. A classification of this 
type was introduced in Refs. £,|13|], for atomic- molecular 
BECs and single-component BECs, respectively, in terms 
of four single component states. These states are the 
OS - on-site symmetric and OA - on-site antisymmetric 
modes, which are respectively symmetric and antisym- 
metric with respect to a minimum of the OL, hereafter 
referred as lattice site, and the IS - inter-site symmet- 
ric and lA - inter-site anti- symmetric modes which are 
respectively symmetric and antisymmetric with respect 
to a maximum of the OL (the center of the symmetry is 
the middle point between two lattice sites). Notice that 
these symmetry properties are the same as for intrinsic 
localized modes of nonlinear lattices, as one could ex- 
pect from the tight-binding approximation based on the 
properties of Wannier functions [14] and reducing the pe- 
riodic GP equations to discrete lattices [7|,|l5|. Except for 
the OS solution, which is always stable, the other types 
of gap solitons are usually unstable in the one compo- 
nent case [13]. However, from ([9]) one can expect that 
gap solitons which have the same symmetry properties 
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in both components should exist and may be stable. The 
respective modes will be characterized by a double in- 
dexing, like for example OS-OS, OS-OA, etc., the first 
and second pairs of letters referring to the first and sec- 
ond components, respectively. The modes OS-OS, OA- 
OA, IS-IS, and lA-IA will be called the " equal- symmetry " 
two- component modes^ while modes OS-OA, IS-IA, etc. 
will be called mixed symmetry modes, with an evident 
meaning in the notation. 



C. Spinor condensates 

Although our results have general character, and can 
be applied to any binary mixture of BECs, either spinor 
or of different species, to be specific we focus on spinor 
BECs. Such condensates are available experimentally. 
As example we mention mixtures of the hyperfine states 
\F = l^rriF = —1) and \F = 2,mi? = 1) as well as of 
the states \F = 2, m^? = 1) and \F = 2, = 2) of 
^^Rb atoms produced in [16] and [17], respectively. For 
a spinor mixture mi = m2 = m, M = m/2, and the 
potential ([TQ|) the equations ([5]) are simplified 

/ii^i = --^ - Focos(2X)$i 

+ (Xi|^i|'+Xl^2|')^i (11a) 

/i2^2 = --^ - FOCOS(2X)$2 

+ (X|^l|'+X2|^2|')^2 (lib) 

where nonlinear coefficients are given by 
Xi = l^~^aii/|aii|, X2 = ^a22/ \an\, and 
X = 2ai2n/[\aii\{n'^ + 1)] with Q = ^/uJ^JuJ^. 

In order to fix parameters we recall that the thermo- 
dynamical stability of a homogeneous binary mixture is 
determined by the condition [18] A > 0, where 

A = XiX2-x'- (12) 

Using this fact as the reference point we will consider the 
following three situations A > 0, A < and A = (we 
however emphasize that, generally speaking, these condi- 
tion should not be considered as the stability /instability 
conditions for localized modes). Before going into details 
of numerical simulations, we observe that interchange 
of the regimes with different A can be experimentally 
achieved either by means of Feshbach resonance affect- 
ing scattering lengths aij or by means of using transverse 
magnetic trap, i.e. geometry of the trap, resulting in dif- 
ferent (jj for different components, say, in ft = V2 for 
the mixtures of hyperfine states \F = 2, mi? = 1) and 
|F = 2, mi? = 2) of ^'^Rb atoms in a magnetic trap ^17j . 



III. LOCALIZED MODES OF REPULSIVE BEC 
MIXTURES WITH UNBALANCED NUMBER OF 
ATOMS (LIMIT iVi > iV2). 

A. Defect modes 

Let us consider the limit A^2 ^ 0. For the first step, 
neglecting contribution of the second component, |<^>2p, 
one can solve equation (|llap in the form 

Mi^r = - K,cos(2X)$f) + Xi|#P*n^l3) 

to find localized modes for a given value of the chemical 
potential jui. In what follows we focus only on symmet- 
ric modes belonging to the lowest branch of the solutions, 
and therefore referred to as the fundamental mode. As 
the second step, in the same approximation, i.e. neglect- 
ing |^2p, one can consider (|llbp with l^i^^P found from 

m- 

^^'''°' = -l'-^ + ^^ff(^)^'^' (14) 
where the effective potential is defined as: 

VeffiX) = -Focos(2X) + x|$^°^p. (15) 

After finding discrete eigenvalues /i2 for a given /ii, one 
can scan the gap for the first component varying /ii and 
obtaining all eigenvalues /i2(Mi)- In Fig. [1] we present the 
results of this procedure for the cases A > and A < 0. 




FIG. 1: Discrete levels in the first lowest gaps /ii,2 G 
[—0.055, 0.92]. The insets show zoom of branches B (OA sym- 
metry) and C (OS symmetry). The amplitude of the OL is 
Vb = 1. (a) and (b) correspond to different sets of nonlinear 
coefficients Xi-X2-X such as in (a) 1:1.584:1.231 (A > 0) and 
in (b) 1:0.942:0.971 (A < 0, these relations are taken from 
p^). The dashed lines in (a) and (b) correspond to the cross 
sections at /ii = 0.86 and /xi = 0.9, respectively. 

The physical meaning of the above procedure consists 
in the approximation of a localized modes ^i, ^2 in the 
limit N2 ^ 0. In this limit, the wavefunction of the 
second component, ^2, appears as a defect mode and, 
strictly speaking, does not describe the atomic distribu- 
tions at a finite N2. Since however, the solution at N2 ^ 



is known, departing from it and using an iterative pro- 
cedure one can construct the solutions for N2 > 0. This 
approach is based on the fohowing important statement 
(whose arguments are presented in Appendix [A|) : // /ii 
and 112 CLS well as Vqi and Vq2 are fixed, then the local- 
ized modes of (0j are isolated. When /ii,2 varies one can 
speak about families of the nonhnear locahzed modes. 

Before performing the search of locahzed modes, let us 
discuss in more details, the results of the above analysis 
of defect modes. In Fig. [T] three branches of the solu- 
tions are shown for different signs of A. The branches 
A correspond to the case where the chemical potentials 
are approximately equal. In order to understand their 
occurrence let us consider a solution of Eqs. (pTj) in the 
form ^2 = OL^i where a is a constant, which without 
loss of generality can be considered real. Such a solution 
exists for /ii = 112 and = {x~ Xi)/{x~ X2) > and 
below will be referred to as a trivial solution. In order 
to obtain the linear limit (p!3|) . (p!4j) one has to consider 
a ^ 0, what is possible only for xi = X- In this last 
case the branches A would coincide with the diagonals of 
the boxes in Fig. [H In the cases at hand xi ^ what 
explains deviations of the lines A from the diagonals. 

The most interesting results, however, are shown in 
the insets of Fig. [U where one can observe the two other 
branches of the solution of the eigenvalue problem ([T4j) 
- branches B and C. These modes appear in the vicinity 
of the opposite gap edges for the first and the second 
components, as this is illustrated in Fig. [2l^a) for the case 
A > 0: the eigenvalues are located near the top of the gap 
for the first component and in the vicinity of the bottom 
of the gap for the second component. For N2 > (but 
small) these eigenvalues give rise to OS and OA modes 
which are characterized by a small difference in energies. 
Their shapes are shown in FigO 

As we have already mentioned these solutions appear 
as defect modes, what is illustrated in Fig. [2] where we 
show the fundamental localized mode of the first compo- 
nent near the upper edge of the gap (panel (b)), corre- 
sponding potential Veff and the respective energy levels, 
which can be interpreted as defect modes (panel (c)). 
The two eigenvalues B and C can be viewed as a split- 
ting of discrete defect level (similar to the level splitting 
in a double- well potential). 

A number of the defect modes depends on the ampli- 
tude of the potential Vq. In all our cases the energy dif- 
ferences between the modes were very small, what makes 
it convenient to speak about appearance of the pairs of 
the modes. Then, one can determine Vn, such that for 
Vn < Vq < Vn-\-i there exist n pairs of the defect modes 
(simultaneously with the fundamental mode). In particu- 
lar, for the cos-like potential and A > (A < 0) we have 
found numerically that Vi ^ 0.43, V2 ^ 1.42 {Vi ^ 0.61, 
^2 ^2.01). 
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FIG. 2: (a) Band structure and localized levels for the values 
of parameters fixed as the case A > in Fig.O the fundamen- 
tal component at /xi = 0.86 (open circle) has OS profile, while 
the second component (filled circles) is given for both OS and 
OA profiles considered at the same values of /J2 = 0; 0.05, 0.1 
(see corresponding profiles in Fig. [3|. (b) The profile of the 
fundamental mode, $1, obtained in the limit A^2 = 0. (c) The 
shape of the respective effective potential ([15]). 



B. The modes with iVi > iV2. 

We have introduced Eqs. (p!3|) . ([T4|) for the sake of ap- 
proximation of localized modes in the limit N2 <C A^i. 
In order to construct the solutions emerging from the 
branches B and C we use an iterative procedure starting 
with the solutions obtained for A^2 ^ 0. Fig. [3] shows 
examples of the obtained results. 

The mode for the second component bifurcates from 
the zero-solution. Its amplitude grows with the detuning 
of the chemical potential fi2 towards the center of the gap 
[see Fig. [3jc)-(f) and Fig. [2] (a)], what is accompanied by 
the increase of the number of particles 7V2. 

In Fig. [4] we present the dependence of the number of 
particles in the second component, N2 on /i2 deviating 
from the corresponding values B and C (see Fig. [T]). We 
notice that in our procedure of search of the modes the 
profile of the fundamental component also changes with 
/i2- The change, however, is very small, it corresponds to 
the change of the number of particles in the first compo- 
nent of order of 10~^, and therefore is practically invisible 
on the scale of Fig. [Sja), (b). 

To study dynamical stability of the obtained modes we 
performed direct numerical integration of equations (|4]) 
with mi = 7722 and uji = UJ2' As initial conditions we 
used perturbed stationary solutions obtained using the 
iterative procedure, the amplitude of perturbation being 
of the order of 5%. 

First, in Fig. [5] we show the case A > 0. One can 
see that the time evolution of the two-component soliton 
is stable. By switching off the coupling parameter, i.e. 
setting X = one obtains the system of two uncoupled 
GP equations. The dominant fundamental component is 
close, in some sense, to an OS mode of the single GP 
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FIG. 3: Localized OS (panels (a)-(d)) and OA (panels (e) 
and (f)) modes for A > (the left column) and for A < 
(the right column) with the same parameters as in Fig. [T] 
In (a) and (b) the profiles of the first component are shown 
for /ii = 0.86 and /xi = 0.9 (indicated by dashed lines in 
Fig. [fl), correspondingly. In [(c), (d)] OS and in [(e), (f)] OA 
solutions for the second component corresponding to branches 
C and B from Fig. [T] are shown. Thick, thin, and dashed lines 
correspond to /i2 = 0.0; 0.05 and 0.1, respectively. The dotted 
line in (a), (b) shows the periodic potential (thin and dashed 
Unes in (a) and (b) are not distinguishable on the scale of the 
figure) . 




FIG. 4: The number of particles of the second component 
vs 112 for the cases A > (panel a) and A < (panel b) 
corresponding to Fig.fS] In (a) /xi = 0.86 with A^i ^ 3.17 and 
in (b) 111 = 0.9 with A^i ^ 3.42. Here solid and dashed lines 
correspond to second component with OS and OA symmetry. 

equation with the periodic potential and that is why it is 
stable. The second component, representing two in-phase 
humps and originally coupled by the effective potential 
due to the first component, becomes unstable at x = 0- 
The two humps in the absence of first component, which 
the played role of the effective barrier between them, start 
to move towards each other and overlap. However due to 
the large velocity after collision two wave packets start 
to move outwards the center. The described behavior 
is observed at the initial stages of the evolution shown 
in Fig. [6] (a), (b). At latter times, reaching the turning 
point they change the direction of motion towards the 
center. After several periods of damped oscillations the 
one-hump OS mode is created. The out-of-phase distri- 
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FIG. 5: Dynamics of density of the first component in (a), 
(c) and of the second component in (b), (d). In [(a), (b)] and 
[(c), (d)] initial profiles are taken from Fig. [3] [(a), (c)] and 
[(a), (e)], correspondingly, with /xi = 0.86 and 112 — 0.1. 

OS - OS OS - OA 




FIG. 6: The same initial conditions as in Fig. [5] with x — ^- 
The real domain of calculation is [-147r; 147r]. 

bution of the second component, shown in Fig. [6] (d) is 
stable representing a single-component stable OA mode. 

OS - OS OS - OA 




FIG. 7: The dynamics of density profiles in (a), (c) of the first 
component with OS symmetry and in (b), (d) of the second 
component with OS and OA symmetries. Initial distributions 
are as in Fig. [3] [(b), (d)] and [(b), (f)] with /ii = 0.9 and /i2 = 
0.1, respectively. 

Next we investigated the case A < 0. Typical results 
are shown in Fig. [71 where one observes instability of 
the both OS-OS (panels (a) and (b)) and OS-OA (panels 
(c) and (d)) modes. It has been checked that the time 
at which the instability is developed increases when 112 
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FIG. 8: Band structure and bound state levels of the two 
component BEG in OL with all attractive Xii — X22 = X12 = 
—0.5, (left panel) and all repulsive xii — X22 = X12 = 0.5, 
(right panel) interactions. For both cases the amplitude of the 
OL is Vb = 1-5 and the number of atoms is A^i = 1.5, A^2 = 
0.98. Open circles refer to the first component while filled 
squares refer to the second component. Notice that all equal 
symmetry modes are degenerated. Also notice the presence 
of a bound state near the edge of the upper band (left panel 
)and a bound state n the second gap (right panel). 




FIG. 9: Gondensate wavefunctions of two-component BEG in 
OL corresponding to the bound state levels in the left panel of 
Fig. [8] (all attractive case). Thick and thin lines refer to first 
and second component, respectively (open circles and filled 
squares in Fig. [5]). To identify the symmetry of the solutions, 
lattice potentials (scaled by a factor 3) have been shown by 
thin dotted line. 



approaches the lowest edge of the gap, tending to infinity 
as /i2 approaches fiB,c and respectively as 7V2 0. 



IV. LOCALIZED MODES WITH BALANCED 
NUMBER OF ATOMS (LIMIT Ni N2). 

Let us now turn to the study of two-component local- 
ized modes with a balanced i.e. with the same or with 
similar number of particles in each component. These 
modes can be of equal symmetry or of mixed symmetry 
type. Some of the equal symmetry modes were also found 
in Ref. [5]. Balanced mixed symmetry modes have not 
been discussed before. In contrast to the equal symmetry 
solutions, which for many aspects are similar to standard 
modes of the single component GPE, the mixed state 
symmetries are more specific for multicomponent BEG 
systems in OL displaying interesting stability properties. 

Thus, while the OA and the OA- OA modes are unsta- 
ble for the one- and two-component system, respectively, 
the coupling of the OA mode with an OS mode gives rise 
to a stable OS-OA mixed state of two-component BEG 
in OL. The same situation occurs with the IS and I A 
modes. In the following we show that these mixed sym- 
metry solutions exist for both attractive and repulsive 
inter and intra specie interactions and together with the 
on site OS-OS solution, they represent stable localized 
excitation of a two-component BEG in OL. 

In Fig. [8] we depict the lowest part of the band struc- 
ture of Eqs. (pT]) obtained for the cases of all attractive 
Xii < X22 < 0, X < (left panel) and all repulsive 
Xii > O7X22 > 0,x > (right panel) interactions. We 
see that the nonlinearity induces in both cases a number 
of bound states in the band gap structure, located in the 
forbidden zone below the first band for the attractive case 




FIG. 10: Gondensate wavefunctions of two-component BEG 
in OL corresponding to the bound state levels in the right 
panel of Fig. [8] (all repulsive case). Thick and thin lines refer 
to first and second component, respectively (open circles and 
filled squares in Fig. [8|. To identify the symmetry of the 
solutions, lattice potentials (scaled by a factor 3) have been 
shown by thin dotted lines. 



or in the gap between the first two bands for the repulsive 
case. Similar bound states can also form inside the upper 
band gaps (see right panel of Fig. [8]). In this figure bound 
state levels are labeled by symmetry type symbols close 
to the pair of the corresponding levels. Thus the OS-OS 
symbol denotes a bound state consisting of an OS state 
in both components with the corresponding chemical po- 
tentials labeled by the open circle and the full square 
close to the symbol. We find that the two component 
bound states of same symmetry type are degenerate in 
energy (chemical potential) both for all attractive and all 
repulsive interactions. In contrast to the states OS- OS 
and OA- OA which are found to be stable also in absence 
of the interspecies interaction (x = 0), the states lA-IA 
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FIG. 11: Time evolution of the density profiles of the IS-IA 
mode of the all attractive case in Fig.[9l in absence (x = 0, left 
panels) and in presence (x = —0.5, right panels) of interaction 
between the two components. 



|$|2 |$|2 |$|2 |$|2 




FIG. 12: The same as in Fig. [H] but for the IS-IA mode of 
the all repulsive case in Fig. [TOl 



is stable only in presence of interaction between the two 
species and the state IS-IS is unstable also in presence of 
the interspecies interaction (in this case the state decays 
a bit faster without interaction). In Figs. [9l and [TOl we 
depict the eigenstates corresponding to the bound states 
shown in Fig. [8l 

In addition to the equal symmetry two-component 
bound states there are also states with mixed symmetry 
such as the OA-OS and the IS-IA bound states. These 
states are very stable and represent new excitations of 
the two-component system since they can exist only in 
presence of the interspecie interaction. In particular we 
find that the mixed symmetry states OA- OS and lA-IS 
exist and are stable under time evolution only due to the 
interspecie interaction. This is shown in Figs. [TTl and[T2l 
for the IS-IA mode of the attractive and repulsive cases 
in Fig. [9] and Fig. [TOl respectively. We see that in both 
cases modes are stable in both components only in pres- 
ence of interaction and the switching off of the interaction 
makes the second component to decay into an OS mode 



2-rn • \ • 1 • < • rz| Z-T-r 




^ °~ '-'^ ^ OS-OS 
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FIG. 13: Band structure and bound state levels of the two 
component BEG in OL with intra and inter species interac- 
tions fixed as xi — —0.5, X2 = 0.5, X = (left panel) and as 
Xi = — 0.5,X2 = 0.5, X = —0.5 (right panel). Other parame- 
ters are Vb = 1-5, A^i = 1.5, A^2 = -98. Open circles refer to 
the first component, filled squares to the second one. Notice 
that in this case the equal symmetry modes are non degen- 
erated. Also notice the presence of localized modes near the 
edge of the upper bands. 




FIG. 14: Gondensate wavefunctions of two-component BEG 
in OL corresponding to the bound state levels in the right 
panel of Fig. 1131 Thick and thin lines refer to first and second 
component, respectively (open circles and filled squares in 
Fig. [13]). To identify the symmetry of the solutions, lattice 
potentials (scaled by a factor 3) have been shown by thin 
dotted lines. 



with emission of matter. 



We have also investigated the case of mixed interac- 
tions xi < O7X2 > 0,X < (see Figs. I13|14p . In con- 
trast with previous cases where bound state levels were 
located in the semi-infinite gap or in the first band gap, 



9 



in the case of mixed interactions the bound states are ah 
below the first band, in spite of the fact that the sec- 
ond component has a repulsive character. In absence 
of the interaction, however, the bound states of the sec- 
ond component are indeed located in the first band gap 
as expected (see left panel of Fig. [13]). By increasing 
the strength of the interspecies interaction these eigen- 
values get shifted below the first band, this being a con- 
sequence of the attraction between the two components 
(for repulsive interspecies interactions the levels would 
be repelled into the gap making the formation of bound 
states more difficult). Also notice that in contrast with 
the all attractive and all repulsive cases, the equal sym- 
metry localized modes (as well as the mixed symmetry 
ones) are non degenerate (the components of the bound 
state have different chemical potentials), this being due 
to the fact that the unperturbed levels belong to differ- 
ent band gaps. By changing the values of x (by means 
of Feshbach resonances) one can indeed change the posi- 
tion of the bound state levels with respect to the band 
structure, this affecting the properties (in particular the 
stability) of the corresponding localized state. Similarly 
properties are found also for other parameter choices and 
for two-component modes states of higher band gaps. 



V. BREATHERS. 

The method of construction of localized modes based 
on continuation from the limit N2 ^ 0, can be used for 
obtaining rather different solutions of the coupled GP 
equations with a periodic potential, when all the nonlin- 
earity coefficients are equal X = Xi = X2- Indeed, in that 
case the system (j5j) takes the form 



1 S^C^o 

■2^-^ocos(2X)$, 



+x(l^ 



1^ 



3-jl 



(16) 



and is invariant under transformation 



^1 



($i + $2)/V2, 



^2 



(^1 - ^2)/V2. (17) 



Let us now assume that one of the component, is 
much larger than the other one, say Ni ^ N2. By us- 
ing the method described in the previous section we can 
find stationary solutions of the coupled GP equations 
in the form $^(X,T) = ^j{X)e-'^dT (notice that the 
chemical potential of the large-amplitude component is 
located near the top of the gap while that of the small- 
amplitude component is near the bottom of the gap, such 
that /ii — /i2 is of the order of the energy gap. Thus, sub- 
stituting the stationary solutions (X, T) into ([T7j) we 
readily obtain new solutions cpj (X, T) characterized by 
the amplitudes of the same order and by two chemical 
potentials, i.e. representing breathers (in the sense of 
the periodically oscillating density of the condensate). 

In Fig. [15] we show the dynamics of such breathing 
solution which appears to be dynamically stable. 




FIG. 15: In [(a),(c)], [(b), (d)] dynamics of the density of the 
and 1^:^11^ are shown correspondingly. Initial profiles for 
are calculated for the parameters x = I5 Vb = 1, /ii =0.88 
and fi2 = 0.1. In [(a), (b)] the OS and [(c), (d)] the OA profile 
of $2 was used. 
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FIG. 16: Gap-breather mode of OA-OS type for the case 
Xi = X2 = X = 0.5, Vo — 1.5, and for number of atoms 
A^i = 1.5, X2 = 0.98. The top panels show ten snapshots of 
the density profiles of the first (first panel) and second (sec- 
ond panel) component of the gap-breather taken at regular 
intervals of time during one period T = 187.86 of oscillation. 
The bottom panel show the amplitude of the left peaks of 
both components as a function of time. The arrows indicate 
that the oscillation of the two components occurs in opposite 
phase. 



Besides two-components breather modes with unbal- 
anced number of atoms, we also find breather modes with 
mixed symmetry which are very stable even for relatively 
large values of the number of atoms in both components. 
The existence of these modes is expected both from pre- 
vious considerations and from the existence and stability 
of the mixed symmetry modes with balanced number of 
atoms discussed in the previous section. 
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In Fig. [16] a gap-breather mode of the OA-OS type is 
depicted for the case of ah repulsive interactions. We see 
that the oscihations of the two components are out of 
phase and are very regular, persisting for very long time 
without any emission of matter, this indicating that the 
mode is indeed very close to an exact breather solution 
of the two-component GPE system. 

VI. CONCLUSIONS 

In the present paper we have presented a variety of lo- 
calized modes, having stationary localized density distri- 
bution, and breathers, characterized by periodically vary- 
ing density, in binary mixtures of Bose-Einstein conden- 
sates. In particular, we have considered modes having ei- 
ther close or very different numbers of atoms (we termed 
them balanced and unbalanced modes) and character- 
ized by the same or different symmetries of the atomic 
distributions in each components. We have suggested al- 
gorithms of systematic finding families of the modes. We 
exploited the fact that unbalanced modes bifurcate from 
the linear defect modes, the defect being understood as a 
localized excitation of one of the components in the ab- 
sence of another component. This allowed us construct 
exact breather modes in a systematic way. 

We have introduced the classification of the modes on 
the basis of four possible distributions of each of the com- 
ponents, what in the general case results in the modes of 
8 different types, which can be designated as BC1-BC2, 
where B stands either for " O" in the case of the distribu- 
tion centered in a local minimum of the periodic potential 
or for "I" in the case of the mode centered in the maxi- 
mum of the potential, and Cj stands for S or A depending 
whether the wavefunction of the the j-th component is 
even (S) or odd (A). 

Although our analysis was concentrated on spinor 
condensates, it allows for straightforward generalization 
for arbitrary two-component system described by cou- 
pled nonlinear Schrodinger equations with periodic coef- 
ficients. 
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APPENDIX A: SOME COMMENTS ON THE 
SHOOTING METHOD 

In this Appendix we briefly outline the implementation 
of the shooting method for finding localized solutions of 
([5]). Let /ii,2 and Vqi, V02 be fixed and /i^ belong to 
gaps of the potentials VoiU{X). In the limit X ^ 00 the 
equations in ([5]) separate. Both of them become linear 
Hill equations and the asymptotics of ^1^2 (-^) as X 00 
are given by the formula {j = 1,2) 

^j{X)^C,e-^^''Hj{X), (Al) 

where 27r-periodic functions Hi^2 and values 71,2 > 
are uniquely defined by the parameters of equations. 
The parameters Ci,2 specify the solutions of ([5j) which 
tend to zero at X ^ 00. If Ci,2 are fixed the solu- 
tion <I>i^2(X; Ci, C2), subject to the conditions (|Aip at 
X ^ 00, can be found by shooting backward starting 
from large enough value of X where this asymptotics is 
valid. Let us consider, for the sake of definiteness, even 
localized solutions of ([5]). They have to satisfy the con- 
ditions 

dx dx 

This is the system of two equations for two variables Ci 
and C2. Generically, the solutions of this system are 
isolated [in the sense that if Ci and C2 is a solution there 
is some neighborhood of (C^ ; C2) which does not contain 
other solutions of ([5])]. Similar arguments hold for odd 
and also for non-symmetric solutions of ([5]). In this way 
one arrives at the conclusion formulated in the text. 
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